Energy gaps and interaction blockade in confined quantum systems 
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Many-body effects in confined quantum systems pose a challenging problem due to the simul- 
taneous presence of particle-particle interactions and spatial inhomogeneity. Here we investigate 
universal properties of strongly confined particles that turn out to be dramatically different from 
what is observed for electrons in atoms and molecules. We show that for a large class of har- 
monically confined systems, including small quantum dots and optically trapped atoms, many-body 
particle addition and removal energies, and energy gaps, can accurately be obtained from single- 
particle eigenvalues. Transport blockade phenomena are related to the derivative discontinuity of 
the exchange-correlation functional. This implies that they occur very generally, with Coulomb 
blockade being a particular realization of a more general phenomenon. In particular, we predict 
van-der-Waals blockade in cold atom gases in traps. 

PACS numbers: 31.15.Ew, 71.10.-w, 73.21.-b, 32.80.Pj 
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Many-body effects in confined systems of interacting 
quantum particles are a recurring theme in many areas of 
science, ranging from nuclei over molecules to nanostruc- 
tured semiconductors More recently, much interest 
turned to the physics of trapped cold atom gases. With 
an initial focus on Bose- Einstein condensation 0] , today 
also the properties of confined fermionic atoms are 
of much concern, and the ability to manipulate trapped 
atoms recently led to the suggestion of atomtronics 
Confinement is often modeled by harmonic-oscillator po- 
tentials, but in spite of decades of research the many- 
body physics in the microscopic as well as the mesoscopic 
regime is still not fully understood. Even quantities as 
fundamental as particle addition and removal energies 
and energy gaps are hard to calculate if effects of con- 
finement and of particle-particle interactions are of com- 
parable magnitude. 

In electronic-structure calculations, addition and re- 
moval energies and gaps are often calculated from 
density- functional theory (DFT) Q. A large body of 
knowledge has been accumulated on how such calcula- 
tions should be done, and when their results are reliable. 
This knowledge, however, is largely based on the behavior 
of electrons in atoms, molecules and solids. Extrapola- 
tion to other systems is fraught with dangers, and may 
lead astray in many ways. 

In this work, we reassess the calculation of these quan- 
tities in confined systems. A key ingredient of our anal- 
ysis are near-exact ground-state energies, obtained from 
diagonalizing the many-body Hamiltonian, which allow 
an unbiased assessment of approximate schemes. As con- 
crete examples, we consider electrons in small quantum 
dots [H and fermionic atoms in optical traps [1| . 

Surprisingly, we find that accurate particle addition 
and removal energies can be obtained from local-density 
single-particle potentials, which is not at all what one 
would expect from experience with atoms, molecules and 



bulk semiconductors. From addition and removal ener- 
gies we calculate energy gaps, and estimate the effect of 
the derivative discontinuity in confined systems. We re- 
late this discontinuity to Coulomb blockade 0], which 
allows us to adopt a more general view on blockade phe- 
nomena than the usual one, leading, in particular, to 
the prediction of van-der-Waals blockade in systems of 
trapped cold atoms. 

To begin with, we define a few basic terms and con- 
cepts. The particle-removal energy of an iV-particle 
quantum system is defined as 

Er{N) E{N - 1) - E{N) = /(iV) = -e^v(iV), (1) 

and the particle-addition energy is given by 

Ea{N) := E{N) - E{N + 1) = A{N) - -e^v+i(iV + 1). 

(2) 

Here, E{N) denotes the many-body ground-state energy 
of the iV-particle system and eN{M) denotes the A^th 
eigenvalue of the Af -particle system. When applied to 
molecules, Er{N) and Ea{N) become the ionization en- 
ergy, /, and electron affinity. A, respectively. These quan- 
tities are defined as differences between ground-state en- 
ergies pertaining to different systems. In DFT, ground- 
state energies are readily calculated from the ground- 
state density, which in turn is obtained from the Kohn- 
Sham (KS) orbitals. Together with these orbitals, one 
also obtains KS eigenvalues. Although these are some- 
times taken as a zero-order approximation to quasi- 
particle energies, most of them have no physical signif- 
icance. The exception to this rule is the energy of the 
highest occripied state, whose negative is the ionization 
energy 0,[1,0|- The electron affinity can also be obtained 
from a highest occupied KS eigenvalue, albeit that of the 
N + 1 particle system. 

In approximate calculations of ionization energies of 
atoms and molecules by means of the local-density ap- 
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proximation (LDA), or by semi-local improvements, it is 
typically found that values obtained from LDA total en- 
ergies agree much better with experiment than do those 
obtained from eigenvalues. The latter can be off by a 
a huge margin [lO|. For affinities, the situation is even 
worse: the LDA potential often does not even bind one 
additional electron, predicting instead the anion to be 
unstable 
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Commonly, the bad performance 
of LDA eigenvalues for predicting ionization energies and 
electron affinities is attributed to the erroneous asymp- 
totic decay of the LDA exchange-correlation {xc) poten- 
tial. For electrons, the exact xc potential (i.e., the one re- 
producing the exact density) decays as 1/r, whereas the 
potential corresponding to local approximations decays 
as the density itself, i.e., exponentially. The exponential 
decay of the LDA potential means that the outermost 
electrons are not bound strongly enough, rendering their 
energies useless for the prediction of electron addition and 
removal energies. Self-interaction corrections [l3, 14 1 are 
believed to be required to cure this problem. 

The difference between the particle-removal and addi- 
tion energies is the fundamental energy gap 

EgiN) := Er{N) - EaiN) (3) 
= E{N + 1) + E{N - 1) - 2E{N) = £;f ^ + A^, (4) 
= eAr+i(A^) - eN{N) + v^„,+{r) - Va,c.-{r). (5) 

Here Wa;c,±(i") = SE^c[n]/Sn{r)\N±n, with 77 ^ 0+ is 
the exchange-correlation (xc) potential calculated at the 
particle-rich and the particle-poor side of integer particle 
number N. E^^ = eN+i{N) - eN{N) is the KS gap, and 

Arr = E„ ~ E^^ is known as the derivative discontinu- 

I — 1 1 — li 

ity, because in DFT it obeys = Vxc,+ —Vxc,- 15lll6l|. 
Axe describes a gap that opens upon addition of a sin- 
gle particle to the system, but disappears in the absence 
of interactions. Because of these characteristics it is, in 
sohds, identified with the Mott gap [H, [3, [13] ■ 

In quantum dots, the energy gap is traditionally also 
decomposed into two contributions. The single-particle 
contribution Ae describes the effects of quantized energy 
levels due to geometry and confinement. The charging 
energy e^/C is due to many-body effects which discon- 
tinuously raise the energy gap upon addition of one more 
electron. This effect is hard to describe quantum mechan- 
ically, and is therefore typically described phenomenolog- 
ically, by a classical capacitance C [l], . 

Both decompositions of the full gap must add up to 
the same value, so that E^^ + Axe = Ae + jC. In 
the phenomenological approach, no general microscopic 
expressions for Ae and C are given. If Ae is calculated 
from eigenvalues of noninteracting particles, subject only 
to the confining potentials, e^/C accounts for all many- 
body effects and is the only iV-dependent contribution. 
DFT suggests the alternative identification of Ae with 
the KS gap, and e^/C with Axe- In this case, both Ae 
and e^/C depend on interaction and particle number. 



TABLE I: Negative of tlie electron-removal energies of the 
Af-particle dot, obtained exactly {E^'^-^^), from LDA total- 
energy differences (i?^^^'"^^), and from LDA eigenvalues 
are also the negative electron-addition en- 
ergies of the N — 1-particle systems. (Here and below AE and 
ev refer to calculations as total-energy differences and from 
eigenvalues, respectively.) All values are in atomic units. 
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0.15 


1.16 


1.16 


1.24 




0.25 


1.68 


1.68 


1.79 




0.35 


2.14 


2.15 


2.28 


6 


0.15 


1.35 


1.33 


1.41 




0.25 


1.92 


1.91 


2.02 




0.35 


2.44 


2.43 


2.56 


7 


0.15 


1.60 


1.55 


1.62 




0.25 


2.30 


2.24 


2.34 




0.35 


2.92 


2.87 


2.99 



because E^^ is obtained from the eigenvalues of an ef- 
fective potential, containing selfconsistent Hartree and xc 
terms. The alternative identification has the advantage 
that E^^ is routinely obtained from DFT codes. The 
capacitance C then describes the beyond-mean-field con- 
tribution to blockade, which need not vanish for ;3> 1 

Common local and semi-local density functionals do 
not have a discontinuity, and their prediction for the 
many-body fundamental gap is the KS gap, which can 
be wrong by a large margin. In principle, however. Axe 
can be estimated from eigenvalues obtained by separate 
calculation of two different systems : combining the 
right-hand sides of Eq. ^ and ^ with Eq. ([3]), and 
comparing the result with Eq. one finds 



Axe = tN+i{N + 1) - eN+l{N), 



(6) 



which allows one to estimate Axe (defined as a contri- 
bution to the many-body gap) even in situations where 
the functional used to generate the eigenvalues has no 
discontinuity. Below, we calculate the different contri- 
butions to the many-body gap, as well as separate dis- 
continuities and electron addition and removal energies, 
for electrons in quantum dots as well as for harmonically 
confined atoms in optical traps. 

First, we turn to quantum dots. Here, we consider 
Coulomb-interacting electrons, which we treat by exact 
diagonalization and, separately, by the local-density ap- 
proximation to DFT, in the two-dimensional parameter- 
ization of Attaccalite et al. f21]. Results for electron- 
removal energies, electron addition energies and energy 
gaps are summarized in Tables III and ITTl 

Several aspects of these results are surprising, and un- 
expected from experience with atoms, molecules and bulk 
semiconductors. First, we note, in Table [H that electron 
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TABLE II: Energy gaps of quantum dots, multiplied by 10 
for legibility. Seem main text for explanation of symbols. 
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0.15 


2.16 


2.29 


0.547 


1.60 


2.15 




0.25 


3.21 


3.23 


0.862 


2.20 


3.06 




0.35 


3.88 


3.99 


1.07 


2.69 


3.76 


5 


0.15 


1.98 


1.70 


0.110 


1.67 


1.78 




0.25 


2.38 


2.30 


0.122 


2.15 


2.28 




0.35 


3.00 


2.82 


0.127 


2.66 


2.78 


6 


0.15 


2.46 


2.16 


0.654 


1.46 


2.12 




0.25 


3.76 


3.30 


1.22 


2.00 


3.22 




0.35 


4.80 


4.39 


1.84 


2.45 


4.29 



addition and removal energies obtained from eigenvalues 
are quite close to the exact data. For larger particle num- 
bers, addition and removal energies obtained from LDA 
eigenvalues can even be better than those from total- 
energy differences. 

This behavior of addition and removal energies ob- 
tained from eigenvalues is dramatically different from 
a huge body of experience accumulated for atoms, 
molecules and solids, where eigenvalue-based ionization 
energies differ widely from experimentjlOl , and affinities 
sometimes cannot be obtained at all This difference 
shows that much of the problems commonly associated 
with the LDA, or with LDA KS eigenvalues, are not re- 
ally due to the LDA itself. Rather, the explanation must 
lie in a physical difference between harmonically confined 
fermions and molecular systems. 

In particular, and in contrast to common physics lore, 
this explanation cannot be related to the asymptotic de- 
cay of the LDA potential: For harmonic confinement of 
Coulomb-interacting particles, the 1/r behavior of the 
exact xc potential is still observed, both in three (2^ 
and two [23[ dimensions, but the LDA single-particle or- 
bitals and density now decay as Gaussians. As Fig. 24 
of Ref. [l] shows, the self-consistent LDA density agrees 
well with the many-body density, even in the asymp- 
totic region. As a consequence, v^^^{r) now also decays 
as a Gaussian, i.e., even faster than exponentially, which 
could be expected to worsen the performance of the LDA 
eigenvalues, instead of improving it. 

The effective potential, however, is much stronger than 
for atoms, because it contains the external confining po- 
tential. In a harmonic potential, the system is com- 
pletely confined, i.e., any number of electrons is bound 
and there are no continuum states. The absence of con- 
tinuum states in harmonic confinement, which is a re- 
alistic feature of real quantum dots at low energies and 
atoms in optical traps, is behind the improved binding of 
the anion-like states, and thus, by means of Eqs. H]) to 
also behind the other improvements noted. The erro- 



neous asymptotics of the LDA potential are not sampled 
by the confined particles. 

It is also noteworthy that standard proofs 0, d, of 
the identification — ejv(-^) = I all explicitly or implic- 
itly assume that the external potential decays to zero 
as |r| oo, or the closely related fact that the single- 
particle orbitals far away from a finite system decay expo- 
nentially with an energy-dependent exponent. Neither is 
true for harmonic confinement, which grows indefinitely 
as |r| oo, and produces single-particle orbitals that de- 
cay as Gaussians with a universal (energy independent) 
exponent. A generalization of the proof to harmonic con- 
finement has been sketched in the appendix of Ref. [23| j 
and is vindicated by our numerical results. 

A second aspect of the data that deserves further in- 
vestigation is the behavior of the gaps in Table |lTl The 
LDA KS gap e^^"^'^^ , obtained from eigenvalues of an 
occupied and an unoccupied orbital, greatly underesti- 
mates the many-body gap E^^^'^^ . This shows that 
the underestimate arises from the lowest unoccupied or- 
bital, not from any occupied orbital. On the other hand, 
energy gaps E^^^ obtained by adding A^c to E^^^-^^ 
(or, equivalently, from differences of eigenvalues pertain- 
ing to the A^-particle and N + 1-particle KS systems) 
are slightly different from LDA total-energy differences, 
j^LDA,AE ^ being sometimes a bit better and sometimes 
a bit worse, when compared to E^'^'^^ . To within this 
small fiuctuation, many-body gaps of quantum dots can 
thus be obtained from just two self-consistent calcula- 
tions of eigenvalues, instead of from three self-consistent 
calculations of total energies. Note that in these calcula- 
tions the derivative discontinuity, estimated from Eq. ([6|). 
makes a significant contribution to the many-body gap. 

Next, we turn to optically confined atom gases. The 
basic formalism of DFT applies both to electrons in quan- 
tum dots, and cold atoms in traps. The functionals, 
however, are different because electrons interact via the 
Coulomb interaction, whereas atoms interact via the van- 
der-Waals force. In dilute systems of cold atoms, this 
is commonly modeled by a suitable contact interaction. 
Again assuming harmonic external confinement, we per- 
form a Hartree calculation, v^c = 0. The data in Table 
mil show that the atom-removal energies obtained from 
eigenvalues and from total-energy differences are in simi- 
lar good agreement for trapped atoms as for electrons in 
quantum dots, indicating that the key is indeed not the 
asymptotic behavior of the xc potential, but confinement. 

For the trapped atoms, estimates of the derivative dis- 
continuity are also represented in Table IIIII Not un- 
expectedly, the discontinuity arising from the van-der- 
Waals interaction is smaller than that arising from the 
much stronger Coulomb interaction, but still makes a sig- 
nificant contribution to the gap, in particular for larger 
values of the interaction parameter g. 

Based on the similarity to the quantum-dot case, we 
predict that a blockade phenomenon will also occur if 
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TABLE III: Negative of the atom-removal energies of A'^ atoms 
in a harmonic trap, interacting with a contact interaction of 
strength g, calculated as total-energy differences and from 
single-particle eigenvalues, all obtained in the Hartree approx- 
imation. These are also the negative atom-addition energies of 
the N — 1-particle systems. is the Hartree single-particle 
gap, and A^c is the derivative discontinuity, estimated from 
single-particle eigenvalues according to Eq. (|6}. At N = 6 
and g = 1.0 the exact (CI) gap is 0.9376, which agrees well 
with the sum E" + A^^ = 0.9181. 
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1.0 


2.4172 


2.4328 


0.8694 


0.0487 




2.0 


2.7748 


2.8009 


0.7776 


0.0902 




5.0 


3.6457 


3.6913 


0.6113 


0.1871 


10 


1.0 


3.4528 


3.4911 


4.3 X 10"'' 


0.0453 




2.0 


3.8537 


3.9218 


7.8 X 10"'' 


0.0804 




5.0 


4.8632 


4.9963 


1.6 X 10"" 


0.1569 



repulsively interacting atoms are channeled one by one 
through an optical trap loaded with a small number of 
particles and coupled to a reservoir 25| , even if the inter- 
action between the particles is of the van-der-Waals type. 
Interaction-driven blockade is not limited to the Coulomb 
interaction. Van-der- Waals blockade is expected to play a 
key role in transport experiments on confined cold atoms 
|25| . and in atomtronic devices of the type proposed in 
Ref. [3. 

In conclusion, an unexpected, but very favorable sce- 
nario emerges from this analysis: Many-body particle 
addition and removal energies of confined systems can 
be reliably estimated from single-body energies, and the 
many-body gap and its Coulomb-blockade contribution 
can be obtained with relative ease and good precision, 
even from LDA. These features are expected to apply 
universally to fully confined systems. 

Derivative discontinuities give rise to blockade phe- 
nomena, which we expect to be ubiquitous, the Mott in- 
sulator, Coulomb blockade and van-der-Waals blockade 
being just three particular realizations of a very general 
phenomenon. 
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